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ABSTRACT 



Context. Knowledge of the mass function in open clusters constitutes one way to constrain the formation of low-mass stars and brown 
dwarfs as does the knowledge of the frequency of multiple systems and the properties of disks. 

Aims. The aim of the project is to determine the shape of the mass function in the low-mass and substellar regimes in the pre-main 
sequence (27 Myr) cluster IC4665, which is located at 350 pc from the Sun. 

Methods. We have cross-matched the near-infrared photometric data from the Eighth Data Release (DR8) of the UKIRT Infrared 
Deep Sky Survey (UKIDSS) Galactic Clusters Survey (GCS) with previous optical data obtained with the Canada-France-Hawaii 
(CFH) wide-field camera to improve the determination of the luminosity and mass functions in the low-mass and substellar regimes. 
Results. The availability of i and z photometry taken with the CFH12K camera on the Canada France Hawaii Telescope added strong 
constraints to the UKIDSS photometric selection in this cluster, which is located in a dense region of our Galaxy. We have derived the 
luminosity and mass functions of the cluster down to / = 18.5 mag, corresponding to masses of ~0.025 M Q at the distance and age 
of IC4665 according to theoretical models. In addition, we have extracted new candidate members down to ~20 Jupiter masses in a 
previously unstudied region of the cluster. 

Conclusions. We have derived the mass function over the 0.6-0.04 M Q mass range and found that it is best represented by a log-normal 
function with a peak at 0.25-0.16 M Q , consistent with the determination in the Pleiades. 



Key words. Techniques: photometric — stars: low-mass, brown dwarfs - 
— galaxy: open clusters and associations: individual (IC4665) 



- stars: luminosity function, mass function — infrared: stars 



1. Introduction 

The initial mass function (IMF), the number of stars f ormed in 
the g alactic disk and in cluster s, was first defined by_ Salpeted 
(fl955h . Later, iMiller & Scalol (fl979h and IScalol (I1986I) con- 
ducted a more accurate study of the IMF. It represents an impor- 
tant link between stellar and galactic evolution and is one of the 
key tools for our understanding of the theories of star formation 



Send offprint requests to: N. Lodieu 
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Research Council of Canada, the Institut National des Sciences de 
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and the University of Hawaii. 
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the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



in the high-mass, low-mass, and substellar regimes. The current 
knowledge drawn from the latest sets of data available in the so- 
lar neighbourhood and in clusters point s towards a power-law or 
log-normal represe ntation of the IMF dKroupal 120021: IChabrieil 
l2003HLarsonll2005l) . 

Over the past decade, deep optical surveys conducted 
in 50-600 Myr old open clusters have uncovered numer- 
ous young pre-main-sequence stars and brown dwarfs. 
New stellar and substellar members have been identi- 
fied in clusters, including the Pl eiades ( Bouvier et alJ 



1 19981: IZapatero Osorio et alJ Il999t 
Moraux et al. 2003; Bihain etalj SoQfl 



Dobb ie et al. 



2002t 

Lodieu et alJ 12007a 
Barrado v Navascues et alJ 12002 

2001, 



or Per dStauffer et al.l 119991: 

i Lodieu et al.l 120051). IC 239 1 dBarrado v Navascues et al.. _ 
1 20041: ISpezzi et all 120091). M 35 dBarrado v Navascues et al 
200 lL Blan co 1 dMoraux et al l 120071). IC4665 ddeWitetal 
2006; Careil e & Jamesl l2010h . the Hyades dDobbie et al 
2002at iBouvier et al.l 120081). and NGC2547 dNavloretal 



2002; Ueffries et al.l 12004 Ueffries & Oliveiral 120051) . Cluster 
membership in these works is generally assessed by 
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proper moti on (lHamblv et alJ Il999t lAdamsetalJ 120011: 
Mora ux et al] 120011: iDeacon & Hamblvl 12004 lLodieu et all 
2007 aJ), near-infrared phot o metry (Zapatero Osorio et al. 1997; 
Barrado v Navascues et al] 120021 lLodieu etaD l2007al) . and 



optical spectroscopy for spectral typing, chromospheric activity, 

and s urface gravi ty (Steele & Jameson 1995; Martin et al. 

119961: lLodieu et all 120051: iBihain et all l2010h . The / - J and 
I - K optical-to-infrared colours have proven their efficiency 
to weed out background giants and field dwarfs as empha- 
sised by optical-infrared su rveys c arried out in the Pleiades 
(Zapatero Oso rio et al.|[l9 97). a Per dBarrado v Navascues et al.l 
20021: lLodieu et al.M2005l) . and cr Orionis dBeiar et al.l 120011: 
Caballeroetal .1120071). 



The Two Micron All-sky Survey (2M ASS:ICutri et al 12 003) 
and Deep Near-Infrared Survey (DENIS; Epc htein et al.lll997l) 
provide JHK photometry for all sources down to K s ~ 14.5 
mag, enabling us to confirm the membership of selected low- 
mass cluster candidates in numerous regions. However, those 
panoramic surveys are too shallow to probe substellar candi- 
dates, except in the nearest an d youngest clusters and star- 
forming regions dTei et al.ll2002l) . 

The UKIRT D eep Infrared Sky Survey (UKIDSS; 
lLawrence et al. 2007 fl, made of five sub-surveys, is designed to 
reach three to four magnitudes deeper than 2MASS and cover 
several thousand of square deg rees at infrared wavelen gths. The 
UKIDSS project is defined in lLawrence et al.l fc 007) and uses 
the Wide Field Camera (WFCAM; ICasali et al.M20071) installed 
on the UK InfraRe d Telescope (UKIRT) and the Mauna Kea 
Observatory (MKO: | Tokunaga et alj|2002l) photometric system 
described in iHewett et al.l (120061) . The pipeline processing is 
described in Irwin et al. (2009, in prepp an d the WFCAM 
Science Archive (WSA) in lHamblv et all d2~008l) . 

The Galactic Clusters Survey (GCS), one of the UKIDSS 
sub-surveys, aims at probing young brown dwarfs in ten star- 
forming regions and open clusters over large areas in five pass- 
bands (ZYJHK) across the 1.0-2.5 micron wavelength range 
with a second epo ch in K. The main goal is to measure the 
form of the IMF (ISalpetejll955l iMiller & Scalolll979t IScalol 
1986) in the substellar regime to tackle important issues includ- 
ing the formation and spatial distribution of low-mass stars and 
brown dwarfs. Early results for Upper Sco rpius, in the Pleiades, 
cr Orio nis, an d in Pr aesepe are pre sented i n lLodieu et al .1 (12006. 
I2007h[ I2008L l20lTl). ILodieuetaLl d2007al) . lLodieu et alT(l2009l 
and (Baker et al. 2010), respectively. 

Among the open clusters targeted by the GCS is IC4665, 
whose central coordinates are RA. = 17 h 40 and dec = 05°. 
This cluster is of special interest because it is among the few 
pre-main-sequence open clusters with an age between 10 and 
50 Myr. It has suffered little dynamical evolution and the mea- 
sured mass function is close to the real IMF because star- 
formation processes have ended. The cluste r is young with 
a lithium age of 27 Myr dManzi et al .1 12008|). consiste nt with 
the nuclear age estima te (36 M yr; Mermilliod| 119811) . activ 



ity probes (<40 Myr; Montes & Martin 1998), and other in 



dependent studies (30-4 Myr; iLvngall 19951: iDambisI Il999t 
iKharchenko et al] |2005a). Several distance estimates are avail- 
able for IC 4665 . lAbt & Chaffed (1 19671) found 324 pc from 
isochrone fitting, iLvngal d 1 9951) reported 430 pc in his cata- 
logue of open clusters and Mermilliod (1981) derived ~350 



1 More details at www.ukidss.org 

2 Extensive details on the data reduction of the UKIDSS images is 
available at http://casu.ast.cam.ac.uk/surveys-projects/wfcam/technical 
see also llrwin et al .1 (2004) 



pc. Independent studies based on Hi pparcos observ ations de- 
rived a mean distance close to 3 50 pc : IDambisI (Jl999) obtained 
370+100 pc. iHoogerwerf et all (1200 ll) reported 385+40 pc, and 
Kharchenko et al. (2005a) inferred 352 pc. In this paper, we will 
adopt an age of 30 Myr and a distance of 350 pc. Previous studies 
of the cluster provided a list of low-mass members based on their 
spectral types, radial veloci ties, and lithium abundances (for a 
summary, see iProsserll 19931) . A recent deep optical (I,z) survey 
covering the central 3.82 square degrees in t he cluster identifie d 
new candidate members down to 0.03 M dde Wit eta i] l2006l) . 
Near-infrared photometry was extracted from 2MASS with addi- 
tional photometry obtained for a subsample of brown dwarf can- 
didates with the infrared camera on the Canada-France-Hawaii 
(CFH) telescope. 
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Fig. 1. Schematic view of the UKIDSS GCS coverage in IC 4665 
is shown with thick solid lines. We included kno wn high 
proba bility photometric candidates (s tar symbols) f rom Prosseil 
d 19931) . the optical (7,z) candidates Ide Wit et all d2006l) con- 
firmed (filled circles) and rejected (plus symbols) by the GCS 
photometry, as well as new photometric candidates identified 
in the GCS outside the optical coverage (open triangles). A 
m ore detailed view o f the optical coverage is shown in Fig. 2 
of Ide Witetal1(r20Q6l) . 



We present ZYJHK photometry for a ~4.3 square degree 
area centred on IC4665 obtained with the WFCAM wide- 
field imager on UKIRT within the framework of the UKIDSS 
Galactic Clusters Survey to investigate the low-mass and sub- 
stellar mass function. This paper is organised as follows. In 
Sect. [2] we describe the optical observations obtained with the 
CFHT in 2002. In Sect, fj we present the GCS survey. In Sect.g] 
we describe the selection of candidates, revise the membership 
of the optically selected candidates on the basis of their new 
infrared colours, and explain the identification of new cluster 
member candidates in the GCS. In Sect.|5]we derive the luminos- 
ity and mass function for IC 4665 in the low-mass and substellar 
regimes. 
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2. The optical survey 

The optical survey was described in detail in lde Wit et al.l d2006). 
We provide here a brief summary to bring our infrared observa- 
tions (Sect. [3} into context. The open cluster IC 4665 was tar- 
geted within the framework of a deep, wide-field, photometric 
survey of open clusters and star-forming regions aimed at study- 
ing the influence of the environment on tthe mass function. The 
observations, conducted in June 2002, surveyed 3.82 square de- 
grees in the cluster down to a completeness limit of I,z = 22 
mag, corresponding to ~0.03 M , assuming a distance of 350 
pc and an age of 50-100 Myr dProsseril 1993b . We adopt the fol- 
lowing notation throughout the paper: Z (in capital letter) refers 
to the WFCAM Z filter with a narrow bandpass (centred at 0.8 
/mi) and z refers to the CFH12K optical z filter. The selection 
of cluster member candidates in the (/ - z,I) colour-magnitude 
diagram yielded a total of 691 low-mass and 94 brown dwarf 
candidates, assuming the age and distance indicated above. The 
2MASS catalogue was queried to provide infrared photome- 
try for candidates brighter than K s = 14.3 mag. Additional in- 
frared (JHK) photometry was obtained for a subset of 101 can- 
didates with the CFHT infrared camera (CFHTIR). The number 
of candidates selected from the optical colour-magnitude dia- 
gram alone and later rejected on the basis of their I — K colours 
is about 35 % at faint magnitudes, showing the efficiency of the 
near-infrared photometry. The contamination by field stars was 
also estimated from two control fields. The derived mass func- 
tion after correction for statistical contamination was best fit by 
a lognormal with a m ean mass of 0.32 M , comparable to re- 
sults in the Pleiades dMartrn et al l 119981; iDobbieet al l l2002bi; 
iTei et al .1 12001 iMoraux et al.ll2Q03t iLodieu et al.l l2007ab. and a 
Per (Barrado v Navascues et al. I l2002h . 



3. The UKIDSS GCS infrared survey 

We extracted all point sources from the UKIDSS GCS DR8 
(September 2010) that are located within a three degree radius 
around the centre of IC4665 using a similar Str ucture Query 
Langu age (SQL) query to our work in the Pleiades ( Lodie uet al.l 
2007a). In summary, we selected only point sources detected in 
ZYJHK, as well as point sources with JHK photometry, but 
undetected in Z and Y. The final catalogue contains 303,978 
sources distributed in the sky as illustrated in Fig. Q] The cor- 
responding (Z - J,Z) colour-magnitude diagram is presented in 
Fig. [2] for one out of ten point source. All sources fainter than 
Z = 11.5, J — 11.0, and K — 9.9 mag were retrieved but our 
analysis will solely focus on candidates fainter than J ~ 13.3 
mag. This limit is set by the original selec tion in the (/ - z,z) 
colour-magnitude diagram (Figure 7 in lde Wit et al.l20 06) where 
only objects fainter than / = 15 mag were considered, corre- 
sponding to a mass of ~0.6 M according to 30 Myr isochrones 
dBaraffe etal.il 19981; IChabrier et al] 120001) shifted to a d = 350 
pc. The GCS survey of IC4665 is complete down to / ~ 18.8 
mag and K ~ 18.1 mag, respectively, adopting a power law fit 
to the histograms of the numbers of point sources as a function 
of magnitude. We are able to detect objects a magnitude fainter, 
which we adopt as our detection limits. 

The central coordinates of the tiles observed by the GCS in 
IC4665 are listed in Table [1] and the coverage is shown in Fig. 
Q] The observations were made in May 2006 (Table[T|i. The cov- 
erage of the survey is homogeneous over the entire cluster and 
complete in the J — 12-19 mag range. 

The GCS provides five-band photometry. In addition, we 
have computed proper motions from the ~6 year baseline be- 




Fig. 2. (Z - J,Z) colour-magnitude diagram for 4.3 square de- 
grees surveyed in IC4665 by the UKIDSS GCS. One out of 10 
sources is plotted. 

Table 1. Log of the observations. 



Tile" 


R.A." 


dec." 


Date 


1 


17 h 44 ra 


+06°02' 


2006-05-10 


2 


17 h 44 m 


+05° 10' 


2006-05-10 


3 


17 h 45 m 


+06°55' 


2006-05-11 


4 


17 h 47 ra 


+06°02' 


2006-05-10 


5 


17 h 47 ra 


+05° 10' 


2006-05-10 


6 


17 h 49 m 


+06° 56' 


2006-05-11 


7 


17 h 51 ra 


+06°02' 


2006-05-11 


8 


17 h 51 m 


+05° 10' 


2006-05-10 



" Each WFCAM tile is 54 by 54 arcmin aside and was observed in all 
five passbands (ZYJHK). Average seeing measured on the images 
and in each individual filter ranged from 0.75 to 1.0 arcsec 

b Central coordinates (in J2000) of each tile from UKIDSS GCS DR8 



tween the 2M ASS (observation s taken between 26 April and 
12 June 2000; ICutri et al.1 120031) and GCS observations (May 
2006) for all sources brighter than J - 15.5 mag. Furthermore, 
we have near-infrared photometry for 313 candidates (only six 
have no Z-band photomery) out of the 393 optical candidates in 
the 7=17.0-22.0 magnitud e range, i.e. 77-78% of all the faint 
optical candidates whereas Ide Wit et al] (120061) presented (only) 
-band photometry for 101 faint candidates (i.e. ~25%). 



4. Cross-matching with optical data 

4.1. Optical-to-infrared colour-magnitude diagrams 

The depth and coverage of the GCS observations allow us to pro- 
vide infrared counterparts in five broad-band filters (ZYJHK) for 
the large majority of optical candidates down to the complete- 
ness limit of the CFH12K survey (~0.03 M ). As mentioned in 
the previous section, we have photometry for 313 out of the 393 
(77-78%) optical candidates with / = 17.0-22.0 mag. 
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Table 2. Number of optical candidates with GCS photometry 
(opt) for different intervals of Z magnitude after the different 
selection steps. 



Range 


opt 


SI" 


S2 6 


S3'' 


S4 d 


S5 e 


RP 


R2 g 


Z= 14-15 


91 


72 


72 


72 


70 


68 


74.7 


94.4 


Z = 15-16 


158 


139 


139 


139 


138 


137 


86.7 


98.6 


Z = 16-17 


185 


152 


152 


152 


151 


150 


81.1 


98.7 


Z= 17-18 


139 


101 


101 


101 


101 


100 


71.9 


99.0 


Z = 18-19 


42 


28 


25 


24 


24 


24 


57.1 


85.7 


Z = 19-20 


15 


9 


9 


7 


7 


7 


46.7 


77.7 


Z = 20-21 


15 


8 


8 


6 


6 


6 


40.0 


75.0 


Z < 17.2 


479 


390 


390 


390 


386 


382 


79.7 


97.9 


Z > 17.2 


175 


120 


117 


112 


112 


111 


63.4 


92.5 


Z = 14-21 


654 


510 


507 


502 


498 


493 


75.4 


96.6 



" Number of optically selected candidates kept after applying the Z- 
J selection (step 1; SI) 

Number of candidates kept after the Z - K selection (step 2; S2) 
c Number of candidates kept after the Y — J selection (step 3; S3) 
d Number of candidates kept after the Y - K selection (step 4; S4) 
e Number of candidates kept after the J - K selection (step 5; S5) 
' The Rl percentage refers to the S5/opt ratio 
8 The R2 percentages refers to the S5/S 1 ratio 



We have plotted various colour-magnitude diagrams using 
optical and infrared magnitudes to assess th e membership of th e 
optically selected candidates published in Ide Wit et al.l (|2006). 
On the left-hand side of Fig. [3] we show the (z - J,z) and the 
(Z - J,Z) colour-magnitude diagrams. The same diagrams are 
displayed for the z - K and Z - K colours on the right-hand side 
of Fig. [3] The majority of optical candidates rejected by the z—J 
and z - K colours are also discarded by the Z - J and Z - K 
colours, showing the efficiency of optical-to-infrared colours in 
removing interlopers (Fig. [3). 

Additional colour-magnitude diagrams employed to reject 
photometric contaminants include the (Y - J,Y) and (J - K,J) 
presented in Fig. [4] We also examined the (H — KJ-H) colour- 
colour diagram to assess the nature of the contaminants. This 
diagram, not shown here, suggests that the contamination in the 
optical selection originates from a mixture of reddened early- 
type dwarfs and field M dwarfs. 

4.2. Revised membership of optically selected candidates 

To verify th e mem bership of the 785 optical candidates from 
Ide Wit et al.l (2006) using our near-infrared photometry, we 
cross-matched that list with the full GCS catalogue with a match- 
ing radius of five arcsec and found 654 common objects. We 
found that 58 out 785 optical candidates lie in the region outside 
of GCS coverage. Therefore, 785-58=727 objects are within the 
GCS coverage, implying that 727-654=73 sources are not re- 
covered by our GCS selection. The GCS detection of these 73 
sources, located at the position of the optical source, did not 
satisfy the criteria set in the SQL query. We found that 19 of 
the 73 objects are missing photometry in (at least) one of the 
JHK passbands, and the remaining sources do not satisfy the 
point sources criteria (Class and ClassStat parameters in at 
least one filter). This problem is inherent to the SQL query and 
represents about 10% of the total number o f candidates, a valu e 
slightly higher than for crOrionis (4.4-7.8%: lLodieu et all2 009). 
We note that 35 out of the 54 (-65%) sources with ZYJHK mag- 
nitudes from GCS DR8 would remain as photometric member 



candidates after applying the colour selections described in the 
next paragraph. 

We generated various colour-magnitude diagrams (Figs.|3]& 
HJ to design selection cuts based on the position of the 654 opti- 
cal candidates with ZYJHK photometry. Looking at these can- 
didates in the several colour-magnitude diagrams, we generated 
a subsample of candidates by extracting only those to the right 
of lines defined in each diagram below: 

• (Z - J,Z) = (0.72,14.0) to (1.15,17.8) 

• (Z-/,Z) = (1.15,17.8) to (2.3,21.5) 

• (Z-tf,Z) = (1.5,14.0) to (2.0,17.5) 

• (Z-K,Z) = (2.0,17.5) to (3.4,21.5) 

• (Y - J,Y) = (0.45,14.0) to (0.55,17.0) 

• (Y - J,Y) = (0.55,17.0) to (1.15,20.5) 

• (Y - K,Y) = (1.3,14.0) to (1.45,17.0) 

• (Y - K,Y) = (1.45,17.0) to (2.2,20.5) 

• (J - K,J) > 0.8 for J = 13-17 mag 

• (J- K,J) = (0.8,17.0) to (1.3,19.5) 

We are left with 493 optical candidates that satisfy those cri- 
teria, and accordingly we confirm their candidacy as photometric 
members on the basis of the five-band infrared photometry. This 
sample will be refered to as "high-probability members" in the 
rest of the paper. These members have a magnitude range that 
spans Z = 14.2-21.1 mag. Their coordinates, photometry and 
proper motions are available electronically and an example is 
shown in Table |A~T1 They are plotted as filled circles in Figs. [3] 
and|4] The other optical candidates, rejected as cluster members, 
are listed in Table IB. H and plotted as plus symbols in Figs. [3]& 
|4] Tables I A. 1 1 and IB. II list the names of each source according 
to the International Astronomical Union (IAU), the coordinates 
(in J2000), the optical {Iz) and infrared (ZYJHK) magnitudes, 
as well as the proper motions in right ascension and declination 
expressed in mas/yr. We emphasise that proper motions for ob- 
jects fainter than / = 15.5 mag are unreliable because they are 
computed using 2MASS as first epoch and the GCS as a second 
epoch. 

We compared the percentage of contaminants as a function 
of magnitude before and after applying the colour cuts above: we 
observe an increase in the level of contamination (from ~20% to 
~50%) with decreasing luminosity in particular below Z = 17 
mag (corresponding t o 0.15 M H ; T abled, in agreement with the 
conclusions drawn bv lde Wit et al.l (120061) from a limited sample 
of TT-band observa tions a nd the study of two control fields (see 
also iMoraux et al. 2003) in this open cluster located at lower 
galac tic latitude and towards a dense region (see e.g. iProsserl 
1993). The most powerful selection criterion clearly is the Z - J 
colour because the other colour cuts remove on average less than 
10% of the candidates (R2 percentage in Tabled. 

We looked further into the level of contamination expected 
in the GCS data. We selected two control fields with galactic 
coordinates similar to IC4665. These are loc ated about 0.2 de- 
grees outside the cluster radius determined bv lKharchenko et all 
(2005b). The two control fields contain 6610 and 7725 sources 
with ZYJHK photometry, respectively, compared to a total of 
7220 sources in the central 0. 16 square degrees in IC 4665. After 
applying the various selection cuts designed above, we are left 
with 89 candidates in the central region of the cluster compared 
with 57 and 50 candidates in the two control fields. The average 
level of contamination is of the order of 55-60% and is compa- 
rable in the bright and faint magnitude ranges. This level of con- 
tamination is very likely an upper limit because the control fields 
are located very close to the cluster, much closer than the control 
fields used in the optical study (0.2 vs 3 degrees). The estimated 
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Fig. 3. Various colour-magnitude diagrams for ~4.3 square degree surveyed in I C 4665 by the UKIDSS Galac tic Clusters Survey . 
Overplotted are the 30 and 100 Myr NextGen (solid line; iBaraffe et all 1 19981) and DUSTY (dashed line; IChabrier etall feOOO) 
isochrones shifted to a distance of 350 pc. Masses are given on the right-hand side of the plot for the 30 Myr isochrone. Overplotted 
are t he lines for photome tric selection as described in Section l4~2l Filled circles represent the optically selected member candidates 
from lde Wit et al.l (12006) confirmed by infrared photometry. Crosses are optically selected candidates rejected as cluster members 
on the basis of their infrared colours. Top left: (zcfhiik - J,Zcfh\2k) diagram. Top right: (zcfhiik - K,zcfh\ik) diagram. Bottom 
left: (Z - J,Z) diagram. Bottom right: (Z - K,Z) diagram. Note that the 30 Myr NextGen and DUSTY isochrones are plotted on the 
top hand side diagrams with the set of filters indicated on the axis. 



contamination is hi gher than the lev els derived for the Pleiades 4.3. New cluster member candidates 
dMoraux et al.l l2001) and Alpha Per dBarrado v Navascues et 

2002]), most likely because of the lower galactic latitude of Considering the low-mass and substellar sequence defined by 
IC4665. the high-probability members, we have identified new candi- 

dates over the full area covered by the GCS. The same colour 



cuts as defined in Sect. I4.2l returned a total of 2236 sources with 
Z = 14.0-21.25 mag. Of the 2236 sources, 864 are outside the 
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Fig. 4. (Y— J,Y) and (J-K,J) c olour-magnitude dia grams for a ~4.3 square de gree surveyed in IC 4665 by the GCS. Overplotted are 
the 30 Myr NextGen (solid line; Baraf fe et alj|1998l) and DUSTY (dashed line; IChabrier et alj2000h isochrones shifted to a distance 
of 350 pc. Masses are gi ven on the right-ha nd side of each plot for the 30 Myr isochrones. Filled circles represent the optically 
selected candidates from Ide Wit e t al. (2006j) whose membership is confirmed by their infrared photometry. Crosses are optically 
selected candidates rejected as cluster members on the basis of their infrared colours. Overplotted are the lines for photometric 
selection as described in Section l4~2l 



optical coverage whereas 1372 are located in the common area 
between the optical and GCS surveys, including the 493 opti- 
cal candidates confirmed as photometric members by the GCS 
(Sect. l4~2i >. Hence, we have 1372-493=879 potential new can- 
didates (643 are brighter than Z = 17.2 mag and 236 fainter). 

What about these new objects? Many lie in the blue part 
of the selection in the (Z - J,Z) diagram, suggesting that the 
large majority are indeed non-members. This trend is also ob- 
served in the (Y - J,Y) diagram and clearly confirmed in the op- 
tical (i - Zcfht,'cfht) diagram. Furthermore, 635 and 30 (=665) 
out of 879 have poor quality flags in the i and z filters, respec- 
tively. Of th ese 214 cand i dates, only 26 satisfy the colour cri- 
teria set by Ide Wit et al.l (120061) : 23 were not selected by the 
optical study because their quality flags were greater than 1, 
another one was excluded from the criteria set for the flux ra- 
dius parameter in SEXtractor, and the remaining two (UGCS 
J174823.92+055254.3 and 174543.44+062225.5) were over- 
looked bv lde Witet al.l (12006). Therefore, if we consider all 879 
(or 879-2) candidates selected from the GCS but unselected by 
the CFH12K survey, the level of contamination of the GCS se- 
lected in IC4665 would be of the order of 64%(879/1372 or 
877/1372), consistent with the 55-60% estimated from the con- 
trol fields ((Sect. EE). 

Furthermore, we have 2236- 1372=864 new candidates from 
the GCS located outside the optical coverage, including 664 
brighter than Z = 17.2 mag and 200 fainter (Table ICTT i. If we 
consider the level of contamination inferred above (55-64%), we 
should expect around 239-299 probable members among 664 
bright candidates and about 72-90 out of 200 at fainter mag- 
nitudes, but these estimates are purely statistical. Nevertheless, 
we can say that the probability of membership decreases with 
increasing distance from the selection line. 



4.4. Radial distributions 

To estimate the spatial distribution of the candidates selected 
from the UKIDSS data alone, we calculate the number of objects 
per square degree located within an annulus R,R + dR centred on 
the cl uster centre a = 17 ft 46 m 04'\ 6 = 5 fl '38 m 53 s (Ide Wit et al.l 
2006). The corresponding radial profile is shown in the left-hand 
side of Fig ure [5] This distribution is then fitted by a King profile 
(Kind ll962l) plus a constant assuming that the density of con- 
taminants n cont is uniform: 



n(x) = k 



1 



1 



Vl + x yT+xt 



+ n c . 



(1) 



where k is a normalisation constant, x - (r/r c ) 2 and x, = 
(r t /r c ) 2 . The distance from the cluster centre is r, r c the core 
radius and r, the tidal radius. 

To find the best fit minimizing^ 2 , we left all these parame- 
ters free, except for the t idal radius that we for ced to be between 
6.4 and 12.1 pc. Indeed. [Piskunov et al.l d2008l) found the cluster 
tidal radius to be r, = 8.7 ± 2.3 pc or 10.7 ± 1 .4 pc whether they 
measured it or calculated it from the semi-major axis. We thus 
found £ = 2924+189 sq.deg. -1 , r c = 12.1±0.4 pc, r, = 12.1+0.3 
pc and n co „, = 251 + 8 per square degree. The total covered 
area being 6.26 square degrees, this corresponds to 1571+50 
contaminants over 2236 candidates, or a contamination level of 
70.3+2.3%, which is consistent with our estimate using control 
fields. 

We performed the same analysis for the sample of candidates 
common to both the GCS and the CFH12K survey. The best fit 
gives k = 497+40 sq.deg." 1 , r c =4.3+0.8 pc, r, = 12.1+0.8 pc 
and n cont = 71+7 per square degree (see right panel in Figure|5]l. 
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Fig. 5. Radial distributions for the full GCS coverage (left) and the common area between the GCS and the CFH12K surveys (right). 



Again using that the common area is 3.475 square degrees, we 
find a contamination level of 50+5%. 



4.5. Proper motions for the brightest member candidates 

The mean proper motion of the IC4665 cluster is estimated 
to be (p a cos a,//a)=(-0.57+0.30 - 7 40+0.36) mas/yr from 3 
members included in the study of Kharchenko et al. (2005b), 
in agreement with the a verage proper motion o f 13 members 
measured by Hipparcos (Ho ogerwerf et al.l l2001 ). For more de- 
tails, we refer the reader to the discussion in Section 4.1 of 
Ide Wit et alj d2006l) . Hence, the cluster's motion is low and not 
easy to separate from field stars and reddened background giants. 
However, astrometric information can be used to reject photo- 
metric candidates with large proper motions. 

We derived proper motions for all point sources identified in 
IC4665 and brighter than J - 15.5 mag using 2MASS and the 
GCS as first and second epoch, respectively, thanks to the ~6 
year baseline. We estimated the root mean square of the cross- 
match to better than 15 mas/yr for / < 15.5 mag. We selected 
potential proper motion members among the photometric candi- 
dates, adopting a 3cr (corresponding to a total proper motion less 
than 45 mas/yr; Fig . [6]) cut, which should optimize the selection 
dLodieuetal.ll201ll) . 

Assuming those parameters, we are left with 278 out the 316 
optical candidates brighter than / = 15.5 mag confirmed by the 
GCS as members, yielding a 10% contamination. Repeating this 
process for the new GCS candidates lying outside the CFH12K 
area, we would select 460 out of the 501 candidates with J < 
15.5 mag, again suggesting a contamination of 10%. A similar 
result is obtained for the sample of the GCS candidates with op- 
tical photometry that were not selected as potential members by 
the optical survey. To summarise, the overall level of contamina- 
tion by proper motion non-members is of the order of 10%. 

4.6. Fainter cluster member candidates 

We attempted to identify fainter and lower mass candidates us- 
ing sources that are Z-band non-detections. We applied the same 
colour criteria in Y — J, Y — K, and J - K colours as presented 
in Sect. 14.21 However, we limited our selection to J = 18.8 mag, 
the completeness of the GCS towards IC4665. At this magni- 
tude, point sources have a signal-to-noise ratio of the order of 
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Fig. 6. Vector point diagram for all optical candidates confirmed 
as probable members by the GCS photometry (black dots) and 
photometric non-members (crosses). Typical error bars on the 
proper motion are 15 mas/yr. 



6-8. The SQL query returned three potential candidates after re- 
moving those clearly detected on the Z images. However, two of 
them are actually detected on the deep CFH12K images (Table 

suggesting that they do not belong to the cluster. These three 
candidates are listed at the top of Table [3] 

We repeated the same procedure for Z and Y non-detections, 
using the J - K colour as sole criterion. The query returned two 
good candidates after removal of objects located at the edge of 
the detector and those clearly detected on the Z and Y images. 
These two candidates are listed at the bottom of Table [3] 

We repeated those selections in the cluster centre and the two 
control fields, but the few objects returned in each area were all 
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Table 3. Faint YJHK (top) and JHK (bottom) candidates iden- 
tified in the UKIDSS GCS. 



R.A. fl 


dec" 


Y" 


r 


H" 


K" 


17:43:18.53 
17:46:27.11 
17:50:28.78 


+05:14:34.9 
+06:18:08.8 
+05:51:16.9 


19.929 
19.919 
19.379 


18.657 
18.615 
18.265 


18.047 
18.187 
17.436 


17.326" 
17.269* 
16.820 c 


17:42:04.01 
17:48:42.73 


+06:15:56.6 
+05:00:57.6 




18.783 
18.607 


18.023 
17.855 


16.745 f 
17.132 c 



Undetected on the deep CFH12K images 



detected in Z for the YJHK candidates and in Z and Y for the 
JHK candidates. 



5. The luminosity and mass functions 

In this section, we derive the cluster luminosity and mass func- 
tions based on the sample of optical photometric candidates con- 
firmed by the GCS photometry. This sample contains the 493 
cluster member candidates discussed in Sect. 14.21 We did not 
correct the luminosity and mass function for the level of con- 
tamination estimated in section |4"31 because this sample repre- 
sents the highest quality sample extracted to date. We expect a 
low level of contamination in this photometric sample, which 
combines optical and near-infrared photometry: about 10% of 
them may be proper motion non-members (see Section POl l and 
another 5-10% spectroscopic non -members accordin g to our 
spectroscopic study of Upper Sco (lLodieu et al.ll201 H) . An in- 
dependent photome tric and spectroscopic study of Blanco 1 by 
iMoraux et al.l d2007l) suggests that the expected number of con- 
taminants is low after combining optical and infrared photomet- 
ric bands. 



5.1. The luminosity function 

This sample of optical photometric candidates with infrared 
photometry contains 493 sources in the J ~ 13. 1-18.6 mag- 
Adopt ing the lithium age of 27 Myr derived by iManzi et all 
(2008!) for IC4665, the stellar/substellar limit is at K = 15.55 
mag (or J — 16.49 mag) adopting a distance of 350 pc and choos- 
ing the 30 My r available in the N extGen (iBaraffe et al.lll998l) 
and DUSTY ( Chab rier et al.ll2000h models to be most consistent 
with the age of the cluster. Thus, our sample contains ~91% of 
stars and 8-9% (40-46 out of 493) substellar candidates. 

Figure |7]displays the luminosity function, i.e. the number of 
stars as a function of / magnitude. The number of objects in- 
creases up to / ~ 15 mag, then decreases down to J ~ 17.5 
mag, and rema ins flat afterwards. Error bars are Gehrels errors 
(lGehrelsll9 86) rather than Poissonian error bars because the for- 
mer represent a better approximation to the true error for small 
numbers. The upper limit is defined as l+( ~\/{dN + 0.75)) and 
the lower limit as y/(dN - 0.25) assuming a confidence level of 
one sigma. 



5.2. The mass function 

To convert magnitudes into masses we used the NextGen and 
DUSTY models at 30 My r from the Lyon group (Baraff e~et al.1 
ll998HChabrier et al.|[20Q0h . Estimated masses at a given magni- 
tude and directly taken from the mass-magnitude relationship 
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Fig. 7. Cluster luminosity function: number of objects as a func- 
tion of J magnitude (TableHJ. Error bars are Gehrels errors. 

Table 4. Numbers for the luminosity and mass functions. 



AMag" 




AMass 1 ' 


dN" 


dN/dlvF 


dN/dlogM" 


13.160- 


43.776 


13.468 


0.660- 


-0.470 


38 


200.4 


257.7 


13.776- 


14.368 


14.072 


0.470- 


-0.310 


88 


550.0 


486.9 


14.368- 


45.010 


14.689 


0.310- 


-0.200 


101 


918.2 


530.7 


15.010- 


45.605 


15.308 


0.200- 


-0.135 


100 


1538.5 


585.8 


15.605- 


46.155 


15.880 


0.135- 


-0.093 


80 


1904.8 


494.3 


16.155- 


16.826 


16.491 


0.093- 


-0.059 


62 


1550.0 


253.9 


16.826- 


47.296 


17.061 


0.059- 


-0.044 


12 


1333.3 


148.5 


17.296- 


17.690 


17.493 


0.044- 


-0.035 


4 


444.4 


40.2 


17.690- 


48.084 


17.887 


0.035- 


-0.028 


4 


571.4 


41.3 


18.084- 


18.546 


18.315 


0.028- 


-0.022 


3 


500.0 


28.6 



Magnitude ranges in the /-band 
Central magnitudes (J c ) 

Mass ra nges (AMass ) . The mass intervals were chosen identical to 
those in de Wit et al. ( 2006) for direct comparison. We assumed an 
age of 30 Myr and a distance of 350 pc for IC 4665 
Numbers of objects per magnitude interval 

Mass functions in linear (dN/dM) and logarithmic (dN/dlogM) 
scales 



extracted from the models. Samples of masses are shown on 
the right-hand side of the colour-magnitude diagrams in Fig. [3] 
We used those relations to derive masses and construct the mass 
function. The highest mass cluster member candidate has ~0.66 
M Q and the lowest mass object has 0.021 M , the latter being 
outside the magnitude and mass ranges shown in Table|4] We did 
not attempt to correct the mass function for binaries. If IC 4665 
is indeed a pre-main-sequence cluster, the curre nt system mass 
function should reproduce the IM F fairly well (Sal peterlll955l: 
iMiller & Scaldll979HScalolll986h . 

Figure [8] shows the system mass function for IC 4665 over 
the 0.66-0.02 M mass range (filled circles and solid line). The 
mass function in logarithmic scale peaks around 0.25-0.16 M Q 
and decreases down to our completeness limit at ~0.035 M Q 
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Fig. 8 . Cluster ma ss function assuming an age of 30 Myr 
(Manzietal. 2008) and a distance of 350 pc. The conver- 
sion from magnitude s into masses made use of the NextGe n 
and DUSTY models (Baraff e"eTal][T99l IChabrier et al.ll2000h . 
Values are given in Table [4] The vertical dot-dashed lines repre- 
sent the GCS completeness limits in IC4665. 



(dot-dashed line in Fig. [8]). The completeness in terms of mass is 
simil ar to our study in the Pleiades (red open diamonds and solid 
line; [Lodieu et al. 2007a) because the difference in age between 
the Pleiades and IC4665 (125 Myr vs 27 M yr) is partly com- 
pensa ted by the larger distance of IC4665 dHoogerwerf et al.l 
2001). The penultimate point is probably real and lies higher 
than the point before at higher masses. The same trend is ob- 
served in th e previ ous mass function in IC4665 published by 
de Wit et aL ([2006) as well as in the Pleiades mass function. 



Moraux et al. (2007) pointed out that this rise is not caused by 



an intrinsic rise of the mass function but rather reflects a drop 
in the mass-luminosity relationship owing to the onset of dust in 
the atmospheres of low -mass stars (the M7/M8 gap as defined 
bv iDobbie et al] 1200 2c). The last point of our mass function, 
however, suffers from incompleteness and represents a lower 
limit because the mass range is below our completeness limit 
(Table |4j. Similarly, the first point of the mass function is in- 
complete because o f the saturation of the deep CFH12K images 
(Ide Wit et al.ll2006h . 

In Fig. [8] we overplotted the published mass function (blu e 
open triangles and solid line) of IC 4665 bv lde Wit et all (120061) . 
The main difference lies in the choice of the age: we hav e used 
the recent age estimate of ~30 Mvr dManzi et alJl2008|) base d 
on the l ithium method dRebolo et al]ll992t iBasri et al.lll996h . 
whereas Ide Wit et al] d2006h showed the mass function for 5 
and 100 Myr following the earlier work by iProsserl dl993l) . 
We emphasis e the large uncertainty o n the distance of IC 4665 
(385±40 pc; Hoog erwerfet all 1200 lb . This range in distance 
would translate into a shift in mass with small differences in 
the substellar regime (~0.005 M ) but more significant at higher 
masses (several tens of Jupiter masses for masses larger than 0.5 
M G ). We observe a sharper decline of the mass function derived 



from the GCS data than the one derived by|de Wit et a il (120061) . 
suggesting that the level of contamination of the optical survey 
below 0.1 M was higher than t he one estimated b y this study 
despite the correction applied by de Wit et al] d2006l) . 

The mass function derived in our study is very similar 
to the Pleiades mas s function drawn from the GCS dataset 
dLodieu et alJ l2007ah down to 0.05 M Q . Both mass functions 
peak around the same mass (between 0.16 M and 0.25 M ) 
followed by a decrease down to the completeness limit of the 
GCS survey. The difference in mass observed at the M7/M8 
gap (0.05 M Q for the Pleiades vs 0.035 M for IC4665) could 
arise fr om the large uncertain ty on the distance of IC4665 
(40 pc; iHoogerwerf et al. | l200ll) compared to the 5 pc error on 
the Pleiades's distance dJohnsonl 11957b iGatewood et al.l [2000; 
Southwor th et alj|2005l) . Moreover, the Pleiades mass function 
was based on a photometric and proper motion sample, whereas 
the IC4665 dataset is purely photometric because we do not 
have proper motion over the full magnitude range probed by the 
GCS. 



6. Conclusions and future work 

We have presented the outcome of a wide-field near-infrared 
survey of IC4665 conducted by the UKIDSS Galactic Clusters 
Survey. The cross-matching of the GCS survey with a previous 
optical survey to a similar depth leads to a revision of the mem- 
bership of low-mass stars and brown dwarfs in IC 4665 based on 
their optical-to-infrared and infrared colours. The main results 
of this paper can be summarised as follows: 

1. We confirmed a total of 493 photometric candidates iden- 
tified in a pure optical survey as high-probability member 
candidates using five infrared (ZYJHK) filters in addition to 
the two optical (Iz) passbands. 

2. We identified new cluster member candidates in a previously 
unstudied region of the cluster. 

3. We derived the luminosity and mass functions in the J = 
13.1-18.6 mag range, corresponding to masses of 0.66 and 
0.02 M at an age of 30 Myr and a distance of 350 pc. 

4. We found that the mass function is similar in shape to the 
Pleiades. It is best represented by a lognormal function peak- 
ing around 0.25-0. 16 M over the 0.66-0.04 M Q mass range. 

The full area of the cluster will be observed again in K as part 
of the GCS, which will provide additional proper motion mea- 
surements over the full magnitude range to add proper motion 
information to the photometric selection. We are planning an 
optical spectroscopic follow-up with the AAOmega multi-fibre 
spectrogra ph installed on the 3.9-m Australia n Astronomical 
Telescope dLewis et~al ] 120021: ISharp et al.i r2006) to confirm the 
membership of all optical and GCS candidates. This spectro- 
scopic follow-up will provide (i) a full spectroscopic mass func- 
tion from 0.65 to 0.04 M , (ii) a better estimate of the level of 
contamination in the optical and GCS surveys, and (iii) a full 
spectroscopic sequence of 30 Myr-old late-M dwarfs. Moreover, 
optical and near-infrared spectroscopy of IC 4665 members will 
be extremely valuable to bridge t he gap between star-forming 
regions like Upper Sco (5 Myr; iPreibisch & Zinnecker 2002) 
and o lder open clusters like the Pleiades ("125 Mvn lStaufFer et al.l 
1998) to define a temperature scale at 30 Myr. For this reason, 
IC 4665 can serve as a benchmark to study the evolution of the 
binary properties in the substellar regime, investigate the role 
of gravity as a function of age, the evolution of disks, and test 
current evolutionary models. Finally, deeper surveys of IC4665 
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may reveal late-L and T dwarf photometric candidates of inter- 
mediate age that would be ideal targets for photometric and spec- 
troscopic studies with upcoming facilities such as the E-ELT and 
James Webb Space Telescope. 
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Appendix A: Optical candidates confirmed by 
infrared photometry 

Appendix B: Optical candidates rejected by 
infrared photometry 

Appendix C: New photometric candidates identified 
in the GCS 
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Table A.l. Coordinates, optical and near-infrared photometry, proper motions (in mas/yr), luminosities (in L/Lo), and mass (in MO) 
for the 493 optically selected cluster member candidates in IC4665 confirmed as such by their infrared colours extracted from the 
Eighth Data Release of the UKIDSS Galactic Clusters Survey 
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fi a cosS e 
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mag 


mag 


mag 


mag 


mag 


mag 


mag 


mas/yr 


mas/yr 






17:42:05.93+05:24:13.9 


17:42:05.93 +05:24:13.9 


15.579 


15.029 


15.113 


14.535 


13.849 


13.200 


12.831 


-1.80 


-19.56 


-1.29 


0.448 


17:50:46.71+06:21:09.8 


17:50:46.71 +06:21:09.8 


18.167 


17.580 


17.446 


16.912 


16.313 


15.699 


15.336 


-43.42 


25.51 


-2.34 


0.084 



8 Name following the IAU and UKIDSS nomenclatures (see http://www.ukidss.org/archive/archive.html for more details). 

b Coordinates in J2000 from UKIDSS GCS DR8 

1 Optical (I,z) photometry from the CFH12K survey ide Wit et alll2006h 

d Near-infrared (ZYJHK) photometry from UKIDSS GCS DR8 

e Proper motions (mas/yr) in right ascension and declination from the GCS vs 2MASS cross-match only for sources brighter than J = 15.5 mag 



Table B.l. Coordinates, optical and near-infrared photometry, and proper motions for 161 optically selected candidates in IC4665 
classified as photometric non-members based on their infrared colours extracted from the Eighth Data Release of the UKIDSS 
Galactic Clusters Survey 
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+05:52:53.4 
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1147.56 


175040.35+062405.9 


17:50:40.35 


+06:24:05.9 


18.761 


18.066 
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18.525 


17.765 


17.621 


-54.98 


1110.12 



" Name following the IAU and UKIDSS nomenclatures (see http://www.ukidss.org/archive/archive.html for more details). 

>> Coordinates in J2000 from UKIDSS GCS DR8 

' Optical (I,z) photometry from the CFH12K survey jde Wit et al1l2006h 

d Near-infrared {ZYJHK) photometry from UKIDSS GCS DR8 

e Proper motions (mas/yr) in right ascension and declination from the GCS vs 2MASS cross-match only for sources brighter than J = 15.5 mag 



Table C.l. Coordinates and near-infrared photometry for the new cluster member candidates in IC 4665 extracted from the UKIDSS 
GCS DR8 
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14.048 -22.14 17.17 



" Name following the IAU and UKIDSS nomenclatures (see http://www.ukidss.org/archive/archive.html for more details). 
* Coordinates and photometry from UKIDSS GCS DR8 

c Proper motions (mas/yr) in right ascension and declination from the GCS vs 2MASS cross-match only for sources brighter than J = 15.5 mag 



